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Abstract 

We address the question of how the celebrated universality of local correlations for the real 
eigenvalues of Hermitian random matrices of size N x N can be extended to complex eigenvalues 
in the case of random matrices without symmetry. Depending on the location in the spectrum, 
particular large- N limits (the so-called weakly non-Hermitian limits) lead to one-parameter defor- 
mations of the Airy, sine and Bessel kernels into the complex plane. This makes their universality 
highly suggestive for all symmetry classes. We compare all the known limiting real kernels and 
their deformations into the complex plane for all three Dyson indices ft = 1,2,4, corresponding to 
real, complex and quaternion real matrix elements. This includes new results for Airy kernels in 
the complex plane for j3 — 1,4. For the Gaussian ensembles of elliptic Ginibre and non-Hermitian 
Wishart matrices we give all kernels for finite TV, built from orthogonal and skew-orthogonal poly- 
nomials in the complex plane. Finally we comment on how much is known to date regarding the 
universality of these kernels in the complex plane, and discuss some open problems. 



1 Introduction 



The topic of universality in Hermitian Random Matrix Theory (RMT) has attracted a lot of attention 
in the mathematics community recently, particularly in the context of matrices with elements that 
are independent random variables, as reviewed in pQ. The question that one tries to answer is this: 
Under what conditions are the statistics of eigenvalues of N x N matrices with independent Gaussian 
variables the same (for large matrices) as for more general RMT where matrix elements may become 
coupled? This has been answered under very general assumptions, and we refer to some recent reviews 
on invariant [21 [3] and non-invariant pQ ensembles. 

In this short note we would like to advocate the idea that non-Hermitian RMT with eigenvalues 
in the complex plane also warrants the investigation of universality. Apart from the interest in its 
own right, these models have important applications in physics and other sciences (see e.g. [1]). We 
will focus here on RMT that is close to Hermitian, a regime which is particularly important for 
applications in quantum chaotic scattering (see [5] for a review) and Quantum Chromodynamics 
(QCD), for example. In the latter case, the non-Hermiticity may arise from describing the effect of 
quark chemical potential (as reviewed in [BJ), or from finite lattice spacing effects of the Wilson-Dirac 
operator (see [7] as well as [8] for the solution of this non-Hermitian RMT). 

Being a system of N coupled eigenvalues, Hermitian RMT already offers a rich variety of large- N 
limits, where one has to distinguish the bulk and (soft) edge of the spectrum for Wigner-Dyson (WD) 
ensembles, and in addition the origin (hard edge) for Wishart-Laguerre (WL, or chiral) RMT. Not 
surprisingly complex eigenvalues offer even more possibilities. The limit we will investigate is known 
as the weakly non-Hermitian regime; it connects Hermitian and (strongly) non-Hermitian RMT, and 
was first introduced in [9] in the bulk of the spectrum. For strong non-Hermiticity - which includes 
the well-known circular law and the corresponding universality results - we refer to [10] and references 
therein, although the picture here is also far from being complete. 

In the next section we give a brief list of the six non-Hermitian WD and WL ensembles, and indicate 
where they were first solved in the weak limit. There are three principal reasons why we believe that 
universality may hold. First, in some cases two different Gaussian RMT both give the same answers. 
Second, there are heuristic arguments available for i.i.d. matrix elements using supersymmetry [11] . as 
well as for invariant non-Gaussian ensembles using large- ./V factorisation and orthogonal polynomials 
(OP) [12]. Third, the resulting limiting kernels of (skew-) OP look very similar to the corresponding 
kernels of real eigenvalues, being merely one-parameter deformations of them. One of the main goals 
of this paper is to illustrate this fact. For this purpose we give a complete list of all the known Airy, 
sine, and Bessel kernels for real eigenvalues, side-by-side with their deformed kernels in the complex 
plane, where some of our results are new. 

2 Random matrices and their limiting kernels 

In this section we briefly introduce the Gaussian random matrix ensembles that we consider, and give 
a list of the limiting kernels they lead to, for both real and complex eigenvalues. For simplicity we have 
restricted ourselves to Gaussian ensembles in the Hermitian cases, in order to highlight the parallels 
to their non-Hermitian counterparts. 

We begin with the classical WD and Ginibre ensembles in ^2.1( displaying Airy (£ )2.2p and sine 
( §2.3p behaviour at the (soft) edge and in the bulk of the spectrum respectively, as well as their defor- 
mations. We then introduce the WL ensembles and their non-Hermitian counterparts in £ j2.4( in order 
to access the Bessel behaviour ( §2.5p at the origin (or hard edge). The corresponding orthogonal and 
skew-orthogonal Hermite and Laguerre polynomials are given in Appendix [A] and precise statements 
of the limits that lead to the microscopic kernels can be found in Appendix [Bl 
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2.1 Gaussian ensembles with eigenvalues on R and C 

The three classical Gaussian Wigner-Dyson ensembles (the GOE, GUE and GSE) are defined as (13 

O.G/3E 
N 



I dHexp[-PTrH 2 /4] = c Nj p TT I dxjW^Xj) \A N ({x})f . (2.1) 



The random matrix elements Hjy are real, complex, or quaternion real numbers for (3 = 1,2,4 re- 
spectively, with the condition that the N x N matrix H (N is taken to be even for simplicity) is real 
symmetric, complex Hermitian or complex Hermitian and self-dual for /3 = 1, 2, 4. In the first equation 
we integrate over all independent matrix elements denoted by dH. The Gaussian weight completely 
factorises and thus the independent elements are normal random variables; for /3 = 1, for example, the 
real elements are distributed M(0, 1) for off-diagonal elements, and Af(0, y/2) for diagonal elements. 

In the second step above, we diagonalised the matrix H = C/diag(xi, . . . ,xat)J7 _1 where U is an 
orthogonal, unitary or unitary-symplectic matrix for f3 = 1, 2, 4. The integral over the latter factorises 
and leads to the known constants c/v,/3- We obtain a Gaussian weight wp{x) and the Vandermonde 
determinant A^({x}) from the Jacobian of the diagonalisation, 



wp{x) = exp[-/3x 2 /4] , A N ({x}) 



n 



xi 



(2.2) 



KKk<N 



The integrand on the right-hand side of eq. (|2.ip times cn^/Z^ E defines the normalised joint prob 



ability distribution function (jpdf) of all eigenvalues. The fc-point correlation function R 1 ^., which is 
proportional to the jpdf integrated over N — k eigenvalues, can be expressed through a single kernel 
K^f of orthogonal polynomials (OP) for f3 = 2, or through a 2 x 2 matrix-valued kernel involving 
skew-OP for = 1,4: 



R p k 2 (xi,... ,x k ) 



det [Kfr 2 (xi,Xj)] , 
i,j=i,...,k 

^=1,4 



Pf 

i,j=l,...,k 



N 



[Xij Xjj 



K 



-W 



0=1,4 



N 



\X%i Xjj 



(2.3) 



The matrix kernel elements Kn and Wn are not independent of Gat but are related by differentiation 
and integration respectively. These relations will be given later for the limiting kernels. 

The three parameter-dependent Ginibre (i.e. elliptic or Ginibre-Girko) ensembles, denoted by 
GinOE, GinUE, and GinSE, can be written as 



7 Gin^E 



(r) 



dJ exp 



"7/3 



1 



Tr ( J Jt 



dH\ dH2 exp 



7/3 Tr Hi jpTrHl 



l + T 



1 



(2-4) 

with r G [0,1). We use the parametrisation of |lUj . with 7/3=2 = 1 and 7/3=1,4 = \- The matrix 
elements of J are of the same types as for H for all three values of /3, but without any further symmetry 
constraint. Decomposing J = H\ + iH<i into its Hermitian and anti-Hermitian parts, these ensembles 
can be viewed as Gaussian two-matrix models. For r = (maximal non-Hermiticity) the distribution 
for all matrix elements again factorises. In the opposite, i.e. Hermitian, limit, the Ginibre ensembles 
become the Wigner-Dyson cases. The jpdf of complex (and real) eigenvalues can be computed by 
transforming J into the following form, J = U(Z + T)U~ 1 . For (3 = 2 this is the Schur decomposition, 
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with Z =diag(zi , . . . , zn) containing the complex eigenvalues, and T being upper triang ulaj0 

-1 



A' 



^ inUE (r) = J] / d 2 z 3 vJjU( Zj ) \A N ({z})\ 2 , w% 2 (z) = exp 

3=1 JC 



1-T 2 



* 2- 



(2.5) 

For f3 = 1,4 we follow |10| where the two ensembles have been cast into a unifying framework. For 
simplicity we choose N to be even. Here the matrix Z can be chosen to be 2 x 2 block diagonal and 
T to be upper block triangular. The calculation of the jpdf reduces to a 2 x 2 calculation, yielding 



N N/2 

^ o/sE (r)=4^n/^n^M( 

7=1 Jc k=i 



Z 2 k-l,Z 2k ) A N ({z}) 



(2.6) 



where we have introduced an anti-symmetric bivariate weight function. For (3 = 1, this is given by 
J=f =1 (z u z 2 ) =w^ =1 {z 1 )w^ =l {z 2 )(2i5 2 {z l - z* 2 )sign(y 1 ) + d 1 (y 1 )5 1 (y 2 )sign(x2-x 1 )y (2.7) 



and for (3 = 4 by 



erfc 



z — z 



V2(l " r 2 ) 



exp 



2(1 + r) 



' (z 2 + z* 2 ) 



Fp =i (zi,z 2 ) = Wp =4 {zi)wp =4 (z 2 ) (zi - z 2 ) 8{z l - z|), (w^ =A (zj\ = w^ =2 {z) . (2i 



For P = 1, it should be noted that the integrand in eq. (|2.6[) is not always positive, and so a symmetri- 



sation must be applied when determining the correlation functions belowifl. For f3 = 4, the parameter 
iV in eq. (12. 6ft should - in our convention - be taken to be the size of the complex-valued matrix that 
is equivalent to the original quaternion real matrix. 

The correlation functions can be written in a similar form as for the real eigenvalues 



R k c 2 ( z i' z k ) = . det [K^ 2 ( Zi ,z*)] , 

i,.7=l,..., k 



R k,C A ( Z ^---i Z k) = . Pf 

i,j=l,...,fe 



A' 
G 



=1,4 



(Zj , Zj ) Gjy £ (Zj , Zj) 



N,C 



=1,4 



(zj , ) Wj^ £1 ' (zj , Zj ) 



(2.9) 



where the elements of the matrix kernels are related through 

^iV,c' ( z i,Zj) = — / d 2 zA^r C ' (.Zj, z)J^g = i Az, Zj) , 

^jvc' ( z *' z i) = / d 2 z d 2 z' F^ =1 ^{z,i, z)K^^' (z, z')J r g =lj4 (z' , z^) — J^L 1)4 (zj, Zj) . 



(2.10) 



The kernels A^ c (z,z') are given explicitly in Appendix lAl 
For /3 = 1, we can write 



Gvc(2l,Z 2 ) = <y 1 (lte)Gj,r reBj (xi,X2) +G^c (Zi,^), 



(2.11) 



The resulting jpdf of complex eigenvalues for normal matrices with T = at j3 = 2 is the same. 
2 It is, however, possible to write the partition function Z^ lnOE as an integral over a true (i.e. positive) jpdf, by, for 
example, appropriately ordering the eigenvalues; however, such a representation is technically more difficult to work 
with. 
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whereas for (3 = 4 eq. (12. 8ft implies the following relations: 



G /3 N £(z 1 ,z 2 ) = (z 2 - z* 2 )w^ =2 {z 2 )K P N ^{z 1 ,z* 2 ) , 

<= 4 (z 1; z 2 ) = -( Zl - z\){z 2 - zi)w% 2 {z 1 )w^ =2 {z 2 )K^{zt,zi) , (2.12) 

(where we dropped the 'contact' term in the final expression). For this reason, for j3 = 4 we will only 
give one of the matrix kernel elements in the following. 

Note that f3 = 1 is special as the eigenvalues of a real asymmetric matrix are either real or come 
in complex conjugate pairs. Therefore we will have to distinguish kernels (and /c-point densities) of 
real, complex or mixed arguments. 

In order to specify the limiting kernels we first need the behaviour of the mean (or macroscopic) 
spectral density. At large N, and for all three values of {3, the (real) eigenvalues in the Hermitian 
cases are predominantly concentrated within the Wigner semi-circle p sc (x) = (2ttN)~ 1 \/4N — x 2 on 
[—2y/~N, 2y/~N], whereas in the non-Hermitian case, the complex eigenvalues lie mostly within an ellipse 
with half-axes of lengths (1 + r)y/N and (1 — r)y/N, with constant density p e \(z) = (Nn(l — r 2 )) -1 . 
Depending on where (and how) we magnify the spectrum locally, we obtain different asymptotic Airy 
or sine kernels for each f3 = 1, 2,4. In the following we will give all of the known real kernels, see e.g. 
|2j for a complete list and references, together with their deformations into the complex plane. For 
the Bessel kernels which will be shown later we need to consider different matrix ensembles, see §2.41 
below. 



2.2 Limiting Airy kernels on R and C 

When appropriately zooming into the "square root" edge of the semi-circle, the three well-known Airy 
kernels (matrix- valued for f3 = 1,4) are obtained for real eigenvalues. For complex eigenvalues we 
have to consider the vicinity of the eigenvalues on a thin ellipse which have the largest real parts, and 
where the weakly non-Hermitian limit introduced in [14] is defined such that 

a = m y/l - t (2.13) 

remains fixed (see Appendix [B] for the precise details of the scaling of the eigenvalues). This leads to 
one-parameter deformations of the Airy kernels in the complex plane. Whilst the result for (3 = 2 is 
already known [lU [15], our results for /3 = 1, 4 stated below are new [16]. 

R-i, x Ai(xi)Ai'(x 2 ) - Ai'(xi)Ai(x 2 ) f°° , . , , , , . 

§ = 2 : Klf {x u x 2 ) = V 1J K 2 > = / dtAi (xi+t Ai {x 2 + t , 2.14 

x\ - x 2 J 



K^(z!,z 2 ) = -^j=e s + 2 / dt e a t Ai(z 1 +t + — )Ai[z 2 + t + 



GyjTT J V 4 / V 4 J 

(2.15) 

In the Hermitian limit a — > we obtain K^~^.(zi, z 2 ) — > y / 5 1 (yi)5 1 (y 2 )Kp j J~ 2 (xi,x 2 ), with the factor 
in front of the integral in eq. (|2.15p projecting the imaginary parts of the eigenvalues to zero. For 
the integral itself - which is obtained from the limit of the sum of the OP on C given in eq. (|A.4|) 
- the deformation in a is very smooth. The same deformed Airy kernel can be obtained from the 
corresponding WL ensemble eq. (|2.29p [15] with kernel eq. (|A.5|h and is thus universal. 
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poo a 

/3 = 4: G P =\x 1 ,x 2 ) = -\K^f{x l} x 2 ) + \A\{x 1 ) / dt Ai(t), K P =\ Xl ,x 2 ) = —G P =\x 1} x 2 ), 
W^ =4 (x 1 ,x 2 ) = - ds G^ 4 (s,x 2 ) 

/"OO f s 

= -\ J ds J dt(Ai(x 2 + t)Ai(x 1 + s) - Ai(x 2 + s)Ai(xi +tfj, (2.16) 



Gilc(^z 2 ) = ( 2 -!7) 

flQQ / 4 4 \ 

xy (is y dt e l ff2 ( s+t ) Ui^ + s + ^-)Ai(zi+t + ^-) -(zi o^)J- 

The integral in eq. (|2.17|) which is also present in the other two kernel elements, see eq. (|2,12|) . clearly 
reduces to that in eq. (|2.16p in the Hermitian limit, whereas the pre- factors provide the appropriate 
Dirac delta functions. When analysing the Hermitian limit in detail, the real kernel elements Gp~ 4 
and A^r 4 follow from a Taylor expansion of Wj^~£, see [17] for a discussion of the analogous Hermitian 
limit of the Bessel kernel. 

poo / poo \ 

P = 1 : Glf{x 1 ,x 2 ) = -J dtM{x 1 +t)M(x 2 + t)- |Ai(a;i) (l-J dtAi(t)j, 
K^T 1 (x 1 ,x 2 ) = -^G^ 1 (x 1 ,x 2 ) , 

f'OO f*X2 /"OO f'OO 

W^r\ Xl ,x 2 ) = - ds G$r 1 (a,x 2 )-% / dt Ai(t) + \ / ds Ai(a) / Ai(t) 

- | sign(xi - s 2 ) , (2.18) 
^.rcai^i.^) = -e- + > / die ff i Ai^ 1 +t + T jAi(x 2 + t + T 

- | e T2+— Ai^xi + —J (^1 — e^a J dt e */ 2 Ai^ + — 



C<S„„te, ») = - 4rt»)(« - * e# + ^^ Jerfc ( !sT) erfc f M 



1, 



/>00 4 4 

y ^ (f H - 1) Ai(zi + 1 + y) Ai (^2 + x) ' 



tf&^i,^) = lsign(y 2 )G^ com { Zl , z* 2 ) , 

=l( Zl , Z2 ) = ( - 2A( Xl ,x 2 ) + B( Xl )B(x 2 ) + B(x 2 ) - B(x\)\8 l (yijS 1 (y 2 ) 



+ 2i [sign(y 2 )G^ Ical (z*, Zl ) - sign( yi )G^ ieal (zl z 2 ) ^ 

- 2i sign( yi )G^ com {zl, z 2 ) 

- 2i5 2 (z 1 - Z2*)sign(yi) - 5 1 (y 1 )5 1 (y 2 )sign(x 2 - x x ) , 
A(xi,x 2 ) = eV +<T {X \ +X2> y ds J dte^ ff2 ( s+ ')Ai^i + s + ^Ai(x 2 + t + 



poo 4 \ 

B( x ) = e^ + ^l dte^ a2t Ai(x + t + ^-). (2.19) 
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Clearly it holds that G^ iCreal (x\,X2) — > G^i (xi,x 2 ) as a — > 0, whereas the complex part vanishes 
in this Hermitian limit Cr^Tc C om( z i' ^2) — 0. We have also explicitly verified the corresponding limits 
for K l l~^(zi,z 2 ) and W^(z 1: z 2 ). 

2.3 Limiting sine kernels on R and C 

For real eigenvalues the sine kernels are obtained by zooming into the bulk of the spectrum, sufficiently 
far away from the edges. The weakly non-Hermitian limit of the complex eigenvalues introduced in 
[U] is taken such that 

a = N 1/2 VT^ (2.20) 

remains finite (see Appendix|B]for further details). In this limit the macroscopic support of the spectral 
density on an ellipse shrinks to the semi-circle distribution on the real axis, whereas microscopically 
we still have correlations of the eigenvalues in the complex plane. 

The list of the known one-parameter deformations of the sine kernels is as follows: 

P = 2: K^\ Xl , X2 ) = Sm(xl ~ X2) = - l\t cos[(xi - x 2 )t] , (2.21) 

7T{Xl — X 2 ) IT Jq 



KkUzuto) = ^ J dt e- 2 ' 2 co8[(*i - z 2 )t] . (2.22) 



The corresponding spectral density of complex eigenvalues was first derived in [9j using super symmetry, 
and the kernel with all correlations functions in pJJ| using OP, see eq. (|A.4j) . In the Hermitian limit 
cr -> 0, we have K^ c {z 1 ,z 2 ) -)■ v / 8 1 (y 1 )5 l (y 2 )K^ 2 (xi,x 2 ). 

In [11] it was shown using supersymmetric techniques that the same result holds for the microscopic 
density of random matrices with i.i.d. matrix elements for j3 = 1,2. Further arguments in favour of 
universality were added in pL2j for the kernel using large- N factorisation and asymptotic OP. The 
universal parameter is the mean macroscopic spectral density p{xq). 



, ^>/8=4/ \ sin[2(xi — x 2 )] T rB=4t \ ® n$=±( 

• ~ 1 : G sin {X U X 2 ) = - 2 ^ Xi _ X2) , Kin ( X 1> X 2)=q x - i G L (^2), 

rxi-xa „ . I f 1 dt 



PX1—X2 1 pi- (If 

Wi= 4 ( xi ,x 2 ) = J dt Gf= 4 (t,0) = —J — sin[2(xi - x 2 )t] , (2.23) 



The corresponding spectral density of complex eigenvalues was derived in [18] using super symmetry, 
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and the kernel with all correlations functions was derived in [19J using skew-OP leading to eq. (|A,9|) . 



(3 = 1: G^ m \x 1 ,x 2 ) = -K? n 2 (xi,x 2 ) , 

d 1 

K^ i ~ 1 (x 1 ,x 2 ) = ■^—G l ^~ 1 (xi,x 2 ) = - dtt sin[(x 2 - x x )t] , 

OX\ TT Jq 

W^ l {x u x 2 ) = [ 1 2 dt G$T%0) + isign(x! - x 2 ) , (2.25) 
J o 

If 1 2 2 

G finC real( X l' X ?) = ~~ dt COs[(xi - X 2 )t] , 

ft Jo 

G £c,cam(*l>*&) = -2isign(y 2 )erfc (^) K^ c {z u z*) , 
1 f 1 

K ^c( z ^^) = - dtte- a2t2 sm[(z 2 - Zl )t] . (2.26) 
ft Jo 

The kernel elements Gg~^(z\, z 2 ) and K^~^(zi, z 2 ) were derived in [20] using skew-OP, c.f. eq. (|A.ll|h 
The same resulting spectral densities of complex and real eigenvalues were derived previously in |21| 
using a sigma-model calculation, which again indicates universality. It can easily be verified here 
that the Hermitian limit <r > of G s - m ^ rea j {xi,x 2 ) is indeed G^ in l (xi,x 2 ), and that Gf in c com (zi, z 2 ) 
vanishes in this limit. 

2.4 Wishart-Laguerre ensembles with eigenvalues on R and C 

In order to be able to access the Bessel kernels for real and complex eigenvalues as well, we briefly 
introduce the Wishart-Laguerre (or chiral) ensembles (L/3E) and their non-Hermitian counterparts 
(CL/3E). We begin with the former which are defined as 

N 

z]f E = / dWexp[-0TrWW*/2] = c x . iM T I" / dxjwUxj) \A N ({x})\ p . (2.27) 
J jj{ Jr+ 

The elements of the rectangular N x (N + v) matrix W are again real, complex, or quaternion real for 
(3 = 1,2,4, without further symmetry constraints. The integration denoted by dW runs over all the 
independent matrix elements. Because we want to access the so-called hard edge of the spectrum we 
will only consider fixed v = O(l) in the following. The distribution of the positive definite eigenvalues 
Xj of WW^ in the Wishart picture (or equivalently the distribution of the singular values of W in the 
Dirac picture used in QCD) is of the same form as eq. (12. ip . but with different weight functions 

w»(x) = x^^+i)- 1 exp[-Px/2] , (2.28) 

that now depend on (3 in a non-trivial way. Consequently the /c-point correlation functions take the 
same form as in eq. (I2.3p . with the corresponding kernels. 

In analogy to the Ginibre ensembles we define a parameter-dependent family of non-Hermitian 
Wishart-Laguerre (also called complex chiral) ensembles as the following two-matrix model 



(2.29) 



Z^ li3E (t) = J dWdV exp -j—^ Tr (WW ] + V ] V - t(WV + 

with W and V* being two rectangular N x (N + v) matrices. Here we follow the notation of [15J. 
This two-matrix model was first introduced and solved for f3 = 1,2,4 in |22l [23l [24] respectively. 
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For r = the jpdf of all the matrix elements again factorises, and in the opposite limit we have 



Z 



CL/3E 
N 



(r) 



7 L/3E 

as T 



1. Here we are seeking the complex (and real) eigenvalues of the product 



matrix WV iW and V are the off-diagonal blocks of the Dirac matrix that we diagonalise) . Its jpdf 
takes the same form as in eqs. (|2.5p and (|2.6p . but with different weight functions that are no longer 
Gaussian 



v. C / \ 

*V =2 ( Z ) 



\z\ exp 



t{z + z* 
1-r 2 



2\z\ 



I 2i/,C f \ 



(2.30) 



For /3 = 4 the anti-symmetric weight function is defined as in eq. ()2.8|) . For /3 = 1 we explicitly specify 
two functions in the anti-symmetric weight function 



Fp'=i(zi,zz) = Wu(zi, z 2 )siga.(yi) 5 2 (z 1 - z 2 ) + ±h u (xi)h u (x 2 )5(yi)5(y 2 )sign(x 2 - xi) 



h v (x) = 2\x\ 2 exp 



TX 



9u(zi,z 2 ) = 2|ziz 2 |2 exp 



dt 



1-r 2 

T(Z\ + Z 2 



a V 1-r 2 



(2.31) 



1-r 2 



T exp 



(4 + zpt 

2\2 



(1-r 



1 

Tt 



Z 2 ~ Zi 



\Vt 



1-r 2 



which are related by g v (z, z*) — > h v (x) 2 as y — > 0. We give the corresponding kernels in Appendix [XI 
In the large- N limit (with v = O(l) fixed), for all three /3 the real positive Wishart eigenvalues 
of WW^ are concentrated on the interval (0, 4/V], with a density p(x) = (2ttN) (4N — x)/x. This 
is a special case of the Marchenko-Pastur density. After mapping to Dirac eigenvalues A = y/x, 



this becomes the same semi-circle distribution as for WD, but with eigenvalues coming in ±A pairs, 
together with v exactly zero eigenvalues making the origin special. The density of the complex Wishart 
eigenvalues has a singularity at the origin; however, after mapping to the Dirac picture, we obtain a 
macroscopic density function that is flat on an ellipse, just as in the Ginibre case. 

We now give a list of all the known Bessel kernels. For real eigenvalues we follow |25[ [26] where 
a most comprehensive list and references can be found. In some cases the parallel between kernels of 
real and complex eigenvalues is more transparent after using some identities for Bessel functions. 



2.5 Limiting Bessel kernels on K and C 

The hard-edge limit is defined by zooming into the origin (see Appendix [B]), where for the complex 
eigenvalues we have to keep a = \J N{1 — r) fixed as in the weakly non-Hermitian bulk limit eq. (|2.20|) . 
The corresponding limiting kernels are given as follows: 

Z(X\ — X 2 ) Jo 

(2.32) 

KL%(^,z 2 ) = ^- 2 K V (M) 2 K v (M) 2 e^ £dtte-^MWzdMW*} ■ 

(2.33) 

It can be shown that K^~ 2 c (zi, z 2 ) \J S 1 (yi)5 x (y 2 )Q(xi)Q(x 2 )K^~ 2 (x\, x 2 ) as a — > 0. The kernel 
of complex eigenvalues was derived in [23]. The same density following from this kernel was obtained 



9 



from a different Gaussian non-Hermitian one-matrix model [27] using replicas, and is in that sense 
universal. 



= 4: K^J(x x ,x 2 ) = 2 I dtt [ dss*[ J 2u +i(2y/xl st)J 2l/+ i(2^ s) - ( Xl x 2 

Jo Jo 

■r/8=4 



G^J?(xi,x 2 ) = -2y/xi J dt J dss 2 ^J 2u (2^/x 1 ~st)J 2u+1 (2^/x 2 ~s) -tJ 2v (2^/x~x~s)J 2v+ i(2*Jx 2 ~st)). 
W^\xi,x 2 ) = 2y/^J dt dss(j 2v {2^st)J 2u (2^ 2 ~s)-(xi^x 2 f), (2.34) 



K L s 4 A z i> z i) 



^ J\t j\sse- 2(72s2(1+t2 ^J 2u (2^r 1 st)J 2u {2^I- 2 s) - ( Zl o 



Z2. 



(2.35) 



The complex kernel was first derived in [24] , whereas the matching in the Hermitian limit - which can 
best be seen when comparing the kernels K^~f c in eq. fl£5g]) with wg- 4 in eq. (|2.34p - is discussed 
in detail in [T7]. One has to take the weight function into account when taking the Hermitian limit. 



G^Jix^x^ 



8y/X!X 2 J 



dss 2 {y/x^ Jv+iis^Jx^J^ST/x^) - (xi <-> x 2 )} 



wL l ( X1 , X2 ) 



dt t 1 ( ^/x{ t ) Jy- 1 ( yfxi t ) 





X2 



dx 



. . — J v {*Jxx) I dsJ v - 2 (s) , 

4 V%i Jj^ 



() G^J(x u x 2 ), 



dsG Bes( s > x z) ~ fSign(xi - x 2 ) 



(2.36) 



G 



Bes,C,com 



r*P =1 f \ 

^Bes.Creal^ 1 ' X2 ) 



1 



2567TCJ 2 



dss 2 e 2(7 s {^/^Jv+i{s^/z{)J v {s^/z 2 )-(zi^z 2 )} , 



(zi,z 2 ) = — 2isign(y 2 )e^ / — e 6*<? 



'dt _^(,2 +42) _A 



4t 



2 V32cr 4 



z 2 



crfcf ^ |y2| W =1 
V 4a 2 / Bes - C 



[sign(x 2 )]^ 



(-*)"/ dy + 



sign^)] 2 Jo 



x 2es^Kz 



\8a 2 J 32^ 



1 _ CT 2 , __ 2o- 
-e MV^i) + jry^+i 



a 



v pi 
■ 2 i ./li 



£\^(crV) J„+i(sv/zT) - a 2 s [E-i-,(a 2 s 2 ) - Ei^(a 2 s 2 ) ) JJsJx? 

2 2 \ 2 2 



(2.37) 



In the final equation above, E n (x) = J dt is the exponential integral. The kernels K^ e ^ c (zx, z 2 ) 

and G%^ c rea i(^i, x 2 ) were derived in [22], whereas the non-commutativity issues of G^~^ c real (si, x 2 ) 
in the weak limit can be found in 128} 29j. There is numerical evidence from studying some examples 
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of non-Gaussian RMT that the density of the real eigenvalues resulting from G^ eB c rea i(a?i, #2) and 
the corresponding distribution of the smallest eigenvalues may be universal |29j . 

The Hermitian limit is much more involved here (compare however K^~^ c (zi, Z2) and K^~J~(zi, Z2)); 

in particular, as N — > 00 and a — > 0, we find that G^~^ c com (zi, Z2) — > and Gg"^ real (xi, X2) — > 

G^~^(xi,X2), and we refer to [28], [29] for more details. 

3 Discussion and open problems 

In this short article we have collected together all the known kernels for RMT with real eigenvalues 
along with all the known and new kernels for RMT with complex eigenvalues at weak non-Hermiticity. 
This comprises the Airy, sine and Bessel kernels of the three Wigner-Dyson and the three Wishart- 
Laguerre ensembles, as well as their non-Hermitian counterparts. In order to highlight the nature of 
this deformation we have used real integral representations for the kernels of real eigenvalues, rather 
than the asymptotic forms resulting from the Christoffel-Darboux identity for (3 = 2 or from the 
rewriting a la Tracy- Widom for f3 = 1,4. The extra exponential factor (and shift for the Airy case) 
in the integral representation of the kernels on C is a very smooth deformation. This makes it very 
plausible that the universality which is very well studied for real eigenvalues extends to the weakly 
non-Hermitian limit for all ensembles, beyond what is already known for (3 = 2. The universality of 
the factor in front of the integral which contains special functions such as the complementary error 
function or modified Bessel function will be more difficult to establish. However, the presence of these 
factors is crucial when taking the Hermitian limit, in projecting the imaginary parts of the eigenvalues 
to zero. 

Whilst we have already mentioned what is known about universality in the weak limit so far, let us 
give some more open problems. To date, a mathematically rigorous derivation of most of the limiting 
kernels on C is lacking, apart from the complex Airy kernel for (3 = 2 [14 . There is no doubt that the 
kernels we have listed and which have been derived using different techniques such as asymptotic OP, 
supersymmetry or replicas are correct. This is based not only on numerical evidence but also, and more 
importantly, on a comparison with complex eigenvalue spectra in physics, see e.g. [6] and references 
therein, where the complex Bessel kernels for (3 = 2 and 4 were successfully compared with complex 
spectra from QCD and QCD-likc theories. Because the latter are field theories and not Gaussian 
RMT this gives a further indication that universality holds in this regime. Preliminary numerical 
investigations with non-Gaussian, non-Hermitian RMT appear to confirm this |29j for j3 = 1. 

A much more challenging problem will be to show the universality of these kernels on C, either by 
going to non-Gaussian potentials of polynomial or harmonic form, or by considering non-Hermitian 
Wigner matrix ensembles, with elements being independent random variables. 

A further reason why we believe that this universality question is important is that some of the 
kernels on C reappear in the same integral form (with real arguments) when looking at symmetry 
transitions between two different Hermitian RMTs, say from one GUE to another GUE, in a corre- 
sponding "weak" limit. Their eigenvalue correlations are also called parametric. For (3 = 2 this fact 
can be observed for the Bessel, sine |30[ [3"T] and Airy [321 [31] kernels. 

Acknowledgements: The organisers and participants of the workshop "Random Matrix Theory 
and its Applications" at MSRI Berkeley, 13th-17th September 2010, are thanked for many inspiring 
talks and discussions. 
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A Finite- N (skew-) orthogonal polynomials and kernels on C 



In this appendix we specify the orthogonal polynomials (OP) and skew-OP as well as their (skew-) 
symmetric scalar products. These may be used to construct the kernels - which we will list for all the 
above matrix ensembles with complex eigenvalues - in terms of which all /c-point correlation functions 
can be expressed, see eqs. (|2.3|) and (|2.9p . We also highlight the relations between expectation values 
of characteristic polynomials on the one hand and (skew-) OP and their kernels on the other, valid on 
both R and C. 

Starting with (3 = 2 we define the monic OP on M (and C) by 




(A.l) 



with squared norms h^, ^ . Because in our examples all moments exist these OP can be constructed 
via the Gram-Schmidt procedure. Alternatively, they can be written as 

p k (z) = (det[z-H]) k = ]F(5ue/ dHdet[z- H}exp[-f3TrH 2 /4] , (A.2) 

and similarly for the GinUE and (C)LUE, replacing the k x k matrix H with J, or the Wishart 
matrices WW^ and WV respectively. In fact, this relation holds for general weight functions. For the 
Gaussian ensembles we obtain Hermite, and for the WL ensembles Laguerre polynomials on M and on 
C. The corresponding kernels are then obtained by summing over the normalised OP (multiplied by 
the weights). Most conveniently, a second relation to characteristic polynomials exists |33| . 

K^(u,v) = w ^ =2 (u)^ =2 (v)i-^-Uet[u- J]det[v- jt]\ , ( A .3) 

"7v-i,c 

which we state here for the GinUE. Correspondingly it holds for the GUE and /3 = 2 WL ensembles, 
and, indeed, for arbitrary weights. We can now give the two (3 = 2 kernels in the complex plane, 
following jJT] and [23] respectively: 



(3 = 2: 




For the skew-OP related to (3 = 1, 4, we have to distinguish between the skew products for complex 
and real eigenvalues. Because the latter are very well known (see e.g. [13]) we will focus on the former, 
which can be written in a unified way [34] 

/ d 2 z 1 d 2 Z2 F^^Azi, z 2 ) det 
Jc 2 ' 

for skew-OP of even-odd degree, and which is vanishing for even-even and odd-odd degree. Here and 
in the following we again choose N even. Once more, the skew-OP satisfying this can be written as 
follows for = 4 [19] and = 1 [M] 

Q P = lA {z) = (det[z- J})^ , Q^f(z) = (det[z- J](z + c + Tr J))^. (A.7) 



P= 



1.4 



(*l) 



o p=1A 

^21+1 
^21+1 i Z 2 



/=1,4 
'fc,(i/)C 



Ski 



(A.6) 
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Note that, for (3 = 4, the matrix J here should be taken as the complex-valued matrix of size 2k x 2k. 
The odd skew-OP in eq. (|A.7|) are defined only up to a constant c times the even skew-OP. The same 
relation holds for real eigenvalues and arbitrary weights, see [34] for references. Moreover, the anti- 
symmetric kernel matrix element K^~^ (sometimes called the pre-kernel) enjoys a similar relation 
to that in eq. (1A.3[) . as was observed for /3 = 4 |17j and (3 = 1 |22j 

K^\u,v) = (u-v)-^ T:r (det[u-J]det[v- J]) „ (A. 



We list the corresponding kernel matrix elements following [19] and [24] respectively 
£-4: %,c^"J- t(1 _ t) ^2 2- 2-(2fc + l)!!(2I)!l 



H^^H^^ < (A - 9) 

^=4 , , _ _ 2 ^ * 2 2*-*-fc!(fc + „)!(2j)! 2fc+2?+1 
^at^cW^- ^(1.^)2 2_y j^(2k + 2v + + i/)! 

x Z*f - (no «)) , (A.10) 

(recalling that N here is the size of the complex-valued matrix that is equivalent to the original 
quaternion real matrix) and [20J and [22] 

1 N ~ 2 r7_L im 

0») = -sebj 5 fr^i t2,+1 MM;) - <» « »)) ■ < A - 12 » 

The other elements of the matrix-valued kernel follow by integration. 

B Large- N limits at weak non-Hermiticity 

In this appendix we will specify the different large- A^ limits that lead to the limiting kernels listed in 
Section [23 Let us emphasise that these are not all of the possible large- A" limits of the above finite- AT 
kernels that one can take. We will give only those limits where (1 — t)N s = a is kept fixed for some 
5 > 0, limits where the degree of non-Hermiticity is weak. The reason is that it is only these particular 
limiting kernels that relate closely to the known universal kernels on R. However, many of the results 
at strong non-Hermiticity (i.e. where r is AMndependent) can be recovered from the weak limit by 
taking a — > oo and rescaling the complex eigenvalues accordingly. 

Soft edge limit: We consider fluctuations around the right end-point of the long half-axis of the 
supporting ellipse, to obtain from eq. (|A.4|) 



i = (l + T)VN + j£r 6 o = nV<VT=F, (B.13) 

K%=l(X 1 +iY 1 ,X 2 + iY 2 )= Y^^K^{ Zl ,z 2 ) . (B.14) 

T-»l 
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The same limit applies to the WL kernel on C given in eq. (1A.5P [15J. By symmetry we expect the 
same limiting behaviour around the left end-point —(1 + r)y/N as well as for u = O(N). The lim- 
iting kernels for f3 = 1,4 are defined in the same way. Note that the eigenvalues in the bulk of the 
spectrum are at strong non-Hermiticity in this limit, since cr s i ne = iV 1//3 cT Airy - > oo as N — > oo. In 
fact, in order to reach weak non-Hermiticity in the bulk we need to consider the following scaling limit. 

Bulk limit: Without loss of generality we consider fluctuations around the origin, being representative 
of the Gaussian ensembles eq. (|2.4j) . On rescaling, we obtain from eq. (|A.4h [9] 

A + i t£t? > " = N^VT^ , (B.15) 



K^(X 1 +iY 1 ,X 2 +iY 2 ) = km ±K^( Zl ,z 2 ) , (B.16) 

r- >1 

and likewise for (3 = 1,4. The macroscopic spectral density collapses onto the real axis, and becomes 
the semi-circle for our Gaussian ensembles. The functions RkC given by the determinant or Pfaffian 
of the rescaled kernels describe the microscopic correlations in the complex plane. 

If we magnify around any other point |xo| < inside the bulk, then we rescale the fluctuations 

z — xq as in eq. (|B. 15|) . The correlations are then universal when measured in units of the local mean 
density 7rp sc (x ). 

Hard edge limit: Whilst we expect that in the bulk of the spectrum the WL and Gaussian ensembles 
(eqs. (|2,29p and (j2.4j) respectively) show the same behaviour, the origin is singled out in the latter 
case. The rescaling here is given by |23] 

z = W +i lN> - = ^ V VT^, (B-17) 
<tc(*i + m,X 3 + iY 2 ) = km — K^%(z u z 2 ) . (B.18) 



n^-od (AN) 



with the Laguerre polynomials in eq. (|A.5|) displaying a Bessel function asymptotic. 

In all three scaling limits the asymptotic kernels are obtained by replacing the sums with integrals 
(the Christoffel-Darboux identity does not hold for OP in the complex plane), and the Hermite and La- 
guerre polynomials by their corresponding Plancherel-Rotach asymptotics (proven for real arguments) 
in the corresponding region. 

An additional problem arises from the integrations with the anti-symmetric weight function J- 
used to obtain the limiting kernel elements G and W. For (3 = 1 these integrals are not absolutely 
convergent, and hence the limit N — > oo and the integration cannot be interchanged. For a detailed 
discussion we refer to 
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